Universal statistics of the local Green function in wave chaotic systems with absorption 
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We establish a general relation between the statistics of the local Green function for systems with chaotic 
wave scattering and uniform energy loss (absorption) and the two-point correlator of its resolvents for the same 
system without absorption. Within the random matrix approach this kind of a fluctuation dissipation relation 
allows us to derive the explicit analytic expression for the joint distribution function of the real and imaginary 
part of the local Green function for all symmetry classes as well as at an arbitrary degree of time-reversal 
symmetry breaking in the system. The outstanding problem of orthogonal symmetry is further reduced to 
simple quadratures. The results can be applied, in particular, to the experimentally accessible impedance and 
reflection in a microwave cavity attached to a single-mode antenna. 
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Statistical fluctuations of physical observables in quantum 
systems with underlying chaotic classical dynamics are the 
subject of a very active field of research in theoretical and ex- 
perimental physics. A considerable progress was underpinned 
by revealing the apparent universality of the fluctuations in 
systems of very diverse microscopic nature, ranging from 
atomic nuclei and Rydberg atoms in strong external fields, to 
complex molecules, quantum dots, and mesoscopic samples 
fjlB]- This universality allows one to exploit the random ma- 
trix theory (RMT) as a powerful tool for a theoretical analysis 
of generic statistical properties of such systems In many 
cases the quantity, which is readily obtained experimentally, 
is the absorption spectrum for transitions from a given initial 
state to highly excited chaotic states at the energy E. Most 
frequently the absorption spectra look practically continuous 
due to both inevitable level broadening and finite experimental 
resolution. Then the relevant statistics are the distribution and 
the correlation functions of the absorption probability |4]. In 
the simplest situation of uniform level broadening V, the prob- 
lem thus amounts to studying statistical properties of the re- 
solvent (Green function) operator G(E) = (E+iT/2-Hy 1 
associated with the random matrix H, which replaces the ac- 
tual chaotic Hamiltonian. In particular, the imaginary part of 
diagonal entries of G(E) is well known in solid state physics 
as the local density of states (LDoSjand in this capacity its 
statistics enjoyed many studies |4j|5l|6j|7l|8|] 

From the experimental point of view the same universality, 
which makes the use of the RMT legitimate, provides one with 
an attractive possibility to employ simple model systems for 
analyzing generic statistics of the fluctuating quantities. Var- 
ious billiards are nowadays proved to be an ideal playground 
for investigating a variety of quantum chaos phenomena; see 
Ref. 1 2] for a general discussion as well as 19t h (111 llll 2lll 311 
for the current problematic. They are usually realized as elec- 
tromagnetic resonators in a form of two-dimensional cavities 
shaped to ensure the chaoticity of internal scattering and cou- 
pled to waveguides (antennas), which are used to inject mi- 
crowaves into the system as well as to collect the output. In 
particular, it turns out that for a single-mode antenna the local 



Green function G(E) (i.e. a diagonal element of G taken at 
the port position) has the direct physical meaning of the elec- 
tric impedance of the cavity which relates linearly voltages 
and currents at the antenna port; see H 1 Oil for a discussion. In 
this way not only the imaginary part but also the real part of 
G turn out to be experimentally accessible quantities whose 
statistics is, therefore, of considerable interest. Inevitable en- 
ergy losses (absorption) must be taken into account properly 
when describing the experiments theoretically 114tll5ll . 

The majority of the experiments is performed in systems 
which are time-reversal invariant (the so-called orthogonal 
symmetry class of the RMT characterized by Dyson's sym- 
metry index (3 — 1). However, non-perturbative analytical 
results are available presently only for systems with no time- 
reversal symmetry (TRS) (j3 = 2 or unitary symmetry class) 
derived by various methods in [H, and for systems with 
spin-orbit scattering ((3 — 4 or symplectic symmetry class) 
|^| . An attempt |6] to provide an expression for the LDoS 
distribution for the (3 = 1 case can not be considered as partic- 
ularly successful, as the final expression was given in a form 
of an intractable fivefold integral. Very recently a heuristic 
interpolation formula at (3 = 1 incorporating exactly the lim- 
iting cases of strong and weak absorption was suggested 001 
to describe reasonably well the available data at moderate ab- 
sorption. Therefore, an exact analytical treatment of the (3=1 
case remains a considerable theoretical challenge. 

In this Letter we present a novel approach to the problem 
which allows us to derive the joint distribution function of the 
real and imaginary parts of the local Green function exactly at 
arbitrary absorption for the whole regime of the gradual TRS 
breaking between the (3=1 and (3 = 2 symmetry classes. 

The Hamiltonian H of the chaotic system gives rise to TV 
levels (eigenfrequencies) characterized locally in the relevant 
range of the energy E by the mean level spacing A. We con- 
sider, as usual, dimensionless quantities expressed in units of 
A, K{E) = (NA/w)G{E), and define the distribution of in- 
terest as follows: 



V(u, v) = (S(u -ReK)5(v + Im K)) 



(1) 
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Angular brackets stand for the ensemble averaging. In such 
units the mean LDoS (v) = 1. The function iK — Z has also 
the meaning of the normalized cavity impedance Z fioll . 

We start with establishing the general relation between the 
joint distribution function Q at finite absorption (assumed to 
be uniform, T > 0) and the energy autocorrelation function 



z_ -iO-K (E-Cl/2-iO) 
1 

z+ + iO-K (E + n/2 + iQ) 



(2) 



of the resolvents of the local Green function Kq at zero ab- 
sorption (r = 0). Distribution Q can be obtained from (0 
by analytic continuation in f2 from a real to purely imaginary 
value ft = iT as follows. Kq(E) is an analytic function of the 
energy in the upper or lower half-plane and can be thus ana- 
lytically continued to the complex values: Kq(E ± iT/2) = 
u^fiv,v > 0. This allows us to continue then analytically the 
correlation function (|2ji from a pair of its real arguments to the 
complex conjugate one: z + = (z_)* = z' + iz", z" > 0. As 
a result, function (|2ji acquires at fl — iT the following form: 



C(z',z") = Cn=ir{z-,z + ) 



1 



oo />oo 

du I dv 



(z' - u) 2 + (z" + v) 2 
V(u,v) 



(V - u) 2 + {z" + vf 



(3) 



The second line here is due to definition Q. To solve this 
equation for V(u, v), we perform first the Fourier transform 
(FT) C(k, z")= J^dz'e^-'Ciz', z") with respect to z' that 
leads to 

r°° - Trp-W(z"+V) 

C(k,z") = dvT(k lV ) — , (4) 



where P(k, v) is the corresponding FT of V(u, v). Being de- 
rived at z" > 0, Eq. (0} can be analytically continued to the 
whole complex z" plane with a cut along negative Re z". Cal- 
culating then the jump of C(k, z") on the discontinuity line 



-V (v > 0), we finally get the following expression 



V(k,v) 



9M 
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[C(k,-v-iO)-C(k,-v + iO)], (5) 



with the Heaviside step function Q(v). The inverse FT of l|5} 
yields the desired distribution. 

This relationship is one of our central results. It resembles 
(and reduces to) the well-known relation between the spectral 
density of states and the imaginary part of the corresponding 
resolvent operator when the case of one real variable is consid- 
ered. In contrast, the case of the distribution function of two 
real variables requires to deal with the two-point correlation 
function. Physically, the latter is a generalized susceptibility 
which describes a response of the system that allows one to 
treat (0 in the sense of a fluctuation dissipation relation: The 
l.h.s. there stands for the distribution (of K) in the presence of 
dissipation / absorption whereas the correlator (of resolvents 
of K) in the r.h.s. accounts for fluctuations in the system, i.e. 
for arbitrary order correlations in the absence of absorption. 

We proceed now with applications. The main advantage 
of the derived relation is that the correlator is a much easier 
object to calculate analytically than the distribution and such 
a calculation for ideal systems at zero absorption has actually 
already been performed in many interesting cases. Let us con- 
sider the chaotic cavity mentioned already. In this case an ex- 
act result for the correlation function (0 has been previously 
obtained by us in Refs. fUEl. Its analytic continuation to 
complex Vl = iT can be represented generally as follows: 



C(z',z") = 



1 



Z' 2 + 0" + l) 2 



d 2 



cf 2 



dz' 2 dz' 1 



Hi)- 



(6) 

where it is important that the function T(x) depends on z' and 
z" only via the scaling variable x = (z 12 + z" 2 + l)/2z" > 
1. Its explicit form depends on the symmetry present (e.g. 
preserved or broken TRS), the following common structure 
being however generic: 



/l />oo t> oo 

dXoj dXij dA 2 /({A})e^ AlA2 - A °)/ 2 

Here, the dimensionless parameter 7 = 2irT/A (i.e. absorption width T in units of the mean level spacing A) accounts for 
the absorption strength. The real function /({A}) is the only symmetry dependent term. In the crossover regime of gradually 
broken TRS it can be represented explicitly as follows JI7II : 

/({A}) = {(1 - Ag)(l + e- 2Y ) (A? - Al)(l - e- 2Y ) + % 2 ft[(l - + X 2 (l e^)]} > W 

I 

with IZ = Aq + Af + A| — 2AoAiA 2 — 1 and Y = y 2 (l — Aq), where y denotes a crossover driving parameter. Physically, 



(x + Xp) 2 

(i + A 1 A 2 ) 2 -(A2-l)(Ai-l) 



(7) 
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y 2 ~ SEy/A is determined by the energy shift SE y of energy 
levels due to a TRS breaking perturbation (e.g., weak external 
magnetic field in the case of quantum dots). Such an effect 
is conventionally modelled within the framework of RMT by 
means of the "Pandey-Mehta" Hamiltonian 11811 . H = Hs + 
i{y/y/N)HA, with Hs (Ha) being a random real symmetric 
(antisymmetric) matrix with independent Gaussian distributed 
entries. The limit y — ► or oo corresponds to fully preserved 
or broken TRS, respectively. 

Now we apply relation (|5} to Eq. (|6} and then perform the 
inverse FT to get V(u, v). Relegating all technical details to a 
more extended publication, we emphasize the most important 
points. The nontrivial contribution to the distribution comes 
from the second ("connected") part of the correlation function 
(|6) whereas the first ("disconnected") one is easily found to 
yield the singular contribution 5(u)5(v — 1). A careful anal- 
ysis shows that due to specific x-dependence given by Eq. |7} 
the above described procedure for the analytic continuation 
of the connected part of C(k, z") is equivalent to continuing 
T(x) analytically and taking the jump at x = —x ± iO, with 



1 



2 c 



> 1. 



The nonzero imaginary part F(x) = Im T(— x) is thus deter- 
mined at given x by the integration region B x = {(Ai, A2) | 
1 < Ai, 2 < 00, (AiA 2 - xf < (XI - 1)(A| - 1)}, where 
the square root in |7} attains pure imaginary values. Tak- 
ing into account the identity d 2 F(x)/du 2 + d 2 F(x)/dv 2 = 



-2 d 
dx 1 



(x — \)4-F(x) valid for x ^ 1, we arrive finally at 



V{u,v) 



1 * j dF{*) 



1 



An 2 v 2 dx 



dx 



2irv 2 



Po(x). (9) 



This distribution is easily checked to be invariant under the 
change iK — > 1/iK, meaning physically that the impedance 
and its inverse have one and the same distribution function. 

Such a form of the distribution is completely generic, as all 
the symmetry specific dependencies were not essential for the 
above discussion. It can be shown 1 8 ] by exploiting the well- 
known relation S = pr|^j| = \fre %e between the scattering 
matrix S and the local Green function (known as K function 
in this context), that the representation given by the second 
equality in Eq. (|9} is a consequence of the two following prop- 
erties at the so-called perfect coupling, k — 1 11911 : (i) the 
uniform distribution of the scattering phase 8 E (0, 2tt); and 
(ii) the statistical independence of 9 and the S matrix modu- 
lus. This establishes also a physical meaning of x by relating 
it to the reflection coefficient r, thus Pq(x) being the normal- 
ized distribution of x = tz^. Remarkably, Eq. (|9j relates the 
distribution of the local Green function in the closed system 
to that of reflection in the perfectly open one. 

Both these properties can be verified using the methods of 
Ref. 1 14] but only in the cases of preserved or completely bro- 
ken TRS. Our approach proves it generally for the crossover 
regime at an arbitrary degree of TRS breaking. Taking into 
account our findings, we can bring the final result to the fol- 



lowing form: 

F(x) = 



dXidX 2 f({X})e-^ XlX2 - x ° y2 



(x - A ) 



[(A?-l)(A|-l)-(A 1 A 2 -x)2]i/2 



(10) 



A (a(x + l)/2) 



(8/2 



D 



-a(x+l)/2 



At arbitrary values of the crossover parameter y the obtained 
result can be treated only numerically. Further analytical 
study is possible in the limiting cases of pure symmetries con- 
sidered below. 

The simplest case of unitary symmetry (J3 — 2) |20] is cor- 
rectly reproduced from (II Oi as y — > 00. We find |21]: 

fl.W - f 

where it is convenient for the subsequent use to introduce the 
absorption parameter a = 7/3/2 scaled with the symmetry 
index (3, a-dependent constants being A = e a ~ 1 and B = 
1 + a — e a , and the normalization constant 7V 2 = 1. 

As to the case of orthogonal symmetry (f3 — 1), no gen- 
eral result was available in the literature. Expression (II It 
(with 7V/3 = a/ (AT ((3/2+1, a) + Be~ a ) and T[v,a) 



(11) 



J^dtt v - 1 e- t ) was suggested in Ref. d (see also M) to be 
an appropriate interpolation formula at f3 = 1. It incorporates 
correctly both known limiting cases of weak or strong absorp- 
tion and a reasonable agreement with available numerical and 
experimental data was found in a broad range of the absorp- 
tion strength. We proceed with providing an exact analytical 
treatment of this case which amounts to investigating (|8} and 
(II Oi at y = 0. Fortunately, further simplifications are possible 
if one considers the integrated probability distribution 

W(x) = - x2 ~_ 1 dF } x) = r dx P (x) , (12) 



2tt 



dx 



which is a positive monotonically decaying function by defi- 
nition. To this end, we note that it is useful to switch to the 
parametrization of Ref. 12211 to carry out the threefold integra- 
tion. The latter turns out to yield then a sum of decoupled 
terms and, after some algebra, we have been able to cast the 
result in the following final form: 

x + 1 r 



W(x) 



4.7T 



f 1 (w)g 2 (w) + f 2 (w)gi(w) 



+hi(w)j 2 (w) + h 2 (w)ji(w) 
with auxiliary functions defined as follows: 

/iW = 



w=(x-l)/2 



(13) 



dt 



9i H 
hi(w) 
h{w) 



dt- 



(1 
1 



w\ e 



- 7 t/2 



1 - e-' 1 + r 



-ft/2 



dt 



~\(l + tf/ 2 



dt- 



1 



-*) 

e -7*/2 



[ 7 +(l 



y/t\t-w\ vT+l 



(14) 
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FIG. 1 : The reflection coefficient distribution in chaotic systems in- 
variant under time-reversal at perfect coupling is plotted at moderate 
values of the absorption strength 7 = 27rr/A, where systematic de- 
viations between the exact result drawn from Eq. 1 1 3i (solid lines) 
and the interpolation expression <1 U at /3 = 1 (dashed lines) are 
most noticeable. Symbols stand for numerics done for 10 3 realiza- 
tions of 500x500 random GOE matrices. 

their counterpart with the index 2 being given by the same 
expression save for the integration region t G [0, w] instead 
of [w, 00). For an illustration of our findings the reflection 
distribution P(r) = ( 1 _ r y Po(j¥p) is shown on Fig. 1. 

Finally, we mention that the general case of arbitrary trans- 
mission T < 1 can be mapped 130 onto that of perfect 
one considered so far. The scattering phase 8 is then no longer 
uniformly distributed and gets statistically correlated with x. 
However, their joint distribution P(x, 9) can be found |8] to 
be again determined by Pq(x) as follows: 

P(x, 6) = ^Po(xg - v 7 ^ 2 - 1)(.9 2 - 1) cos 9) , (15) 

with g = 2/T—l. This equation provides us further with dis- 
tributions of the phase and reflection coefficient which were 
recently studied experimentally |9]. 

In summary, we provided a new approach to statistical 
properties of the local Green function in quantum chaotic 
or disordered absorptive systems of any symmetry class. It 
would be highly interesting to extend current experimental 
studies II ill of the crossover regime to check our findings. 
Although the validity of explicit formulas given above is re- 
stricted to the completely ergodic ("zero-dimensional") case, 
it is actually possible to adopt the above method 1 25] to quasi- 
one (or higher) dimensional situations when Anderson local- 
ization effects play already an important role |7]. 
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